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Percolation clusters are probably the simplest example for scale-invariant structures which either 
are governed by isotropic scaling-laws ( "self-similarity" ) or — as in the case of directed percolation 
- may display anisotropic scaling behavior ( "self-affinity" ) . Taking advantage of the fact that 
both isotropic and directed bond percolation (with one preferred direction) may be mapped onto 
corresponding variants of (Reggeon) field theory, we discuss the crossover between self-similar and 
self-afnne scaling. This has been a long-standing and yet unsolved problem because it is accompa- 
nied by different upper critical dimensions: d\ — Q for isotropic, and dj? = 5 for directed percolation, 
respectively. Using a generalized subtraction scheme we show that this crossover may nevertheless 
be treated consistently within the framework of renormalization group theory. We identify the 
corresponding crossover exponent, and calculate effective exponents for different length scales and 
the pair correlation function to one-loop order. Thus we are able to predict at which characteristic 
anisotropy scale the crossover should occur. The results are subject to direct tests by both computer 
simulations and experiment. We emphasize the broad range of applicability of the proposed method. 

PACS numbers: 0520, 0540, 6460 



I. INTRODUCTION 

In describing a large variety of patterns in nature concepts from fractal geometry [jjj have become of increasing 
importance over the last years. The simplest kind of scale invariance is self-similarity, i.e., invariance with respect 
to homogeneous dilation or contraction. A second kind of structures, which frequently appears in growth models, are 
self-affine clusters, which we define as characterized by anisotropic scaling (see Fig. 1). 

Perhaps the most simple growth model which incorporates both of these fractal structures is percolation In 
ordinary percolation, sites or bonds are filled at random with probability p. The percolation process then proceeds 
along paths connecting occupied nearest neighbors. The clusters formed by nearest-neighbor links are self-similar, 
i.e., they display isotropic scaling. In directed percolation || the links between nearest neighbors have a bias in 
one preferred direction, such that the percolation process advances along this direction only. (Often this direction is 
referred to as the time direction t and therefore directed percolation proceeds in the direction of increasing time.) The 
size of the clusters in the preferred direction is characterized by a length scale different from that in the perpendicular 
direction. Fig. 1 depicts typical clusters emerging from isotropic (a) and directed (b) percolation, respectively. 

If percolation in the positive t direction is merely favored with a certain probability with respect to the negative 
t direction, but propagation "backward" in time is still admitted, the situation will be more complicated. If the 
"anisotropy" g is low, one expects almost isotropic scaling behavior in a large region of the phase diagram. However, 
when the critical region near the percolation threshold p c is approached, self-affine scaling will become apparent. 
Similarly, if p = p c , and the percolation clusters are viewed at a tiny length scale, hardly any deviations from self- 
similarity will be noticeable. On the other hand, if one proceeds to larger and larger scales, anisotropic scaling will 
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become more and more important, until finally the asymptotic limit of directed percolation is reached. 

Our aim is to provide a quantitative description of the ensuing crossover features from self-similar to self-affine 
scaling, when either the length scale is varied, or the threshold p c is approached, for a biased percolation problem. 
More precisely, we are going to compute at which characteristic anisotropy scale g this crossover should occur, and 
we shall find that this crossover point depends on the specific scaling law and correspondent critical exponent under 
consideration. These predictions are, of course, subject to direct tests by computer simulations and/or experiments, 
which would thus be most desirable in order to check our results. 

For the issues we have in mind, the main quantity of interest is the pair correlation function (or connectivity) 
G(r2,i"i), which measures the probability that the sites T2 and ri are connected by some path irrespective of the 
other sites in the lattice. Lines of constant G hence describe the average shape of the percolating structure. In 
directed percolation r = (x, t) the pair correlation function has to be causal, i.e., one has to add the restriction that 
to, > t\ for G to be non-zero. 

Percolation problems may be reformulated in terms of certain field theories. Cardy and Sugar [Q have shown that 
directed bond percolation is in the same universality class as Reggeon field theory, which has been studied intensively 
by particle physicists in the 1970s, below five dimensions. This universality was also confirmed numerically for d = 2 
and d = 3 (see Refs. ^j-Q ) . A corresponding mapping for isotropic site-bond correlation onto a related field-theoretical 
model was performed by Benzoni and Cardy valid for d < 6 dimensions. An extension of these models to the 
intermediate case of anisotropic, but not entirely directed percolation is straightforward, and we shall use the ensuing 
field theory for our investigation of the crossover from self-similar to self-affine scaling. 

The problem of directed percolation has been studied by various theoretical methods, such as high-temperature 
expansion for the calculation of the exponents in d = 2 || , e expansion Q , and Monte Carlo simulations 
main focus of all these investigations was the determination of the independent critical exponents. 

There are many processes in nature which can be described in terms of biased percolation. An example is the 
reaction of polymerization with the production of a giant macromolecule (gelation or vulcanization), if it occurs 
under anisotropic external conditions. If the seed macromolecule is washed by a flow of solution containing monomer 
groups, it is obvious that the probability of connection along and against the flow is different jll|]. In general, directed 
percolation can be understood as a prototypical model for the spreading of some influence, such as transport in a 
strong external field jjJI , crack propagation , epidemics or forest fires with a bias JlJ] . One more example is the 
propagation of excitations in the system of neurons or neuron-like automatons. 

Our paper is organized as follows. In the following section we start with an outline of the mapping of percolation 
problems onto "dynamical" field theories |l5f] . We define the specific model under consideration here and comment 
on the different upper critical dimensions which play a role in the limiting cases of isotropic and directed percolation, 
respectively. In Sec. Ill we shall then describe an appropriate renormalization procedure allowing for a detailed 
analysis of the crossover scenario. It comprises a generalization of Amit and Goldschmidt's procedure designed for 
bicritical points fl(f | enabling us to treat the fixed points with different upper critical dimension within a unified 
renormalization scheme. The fourth section is devoted to the solution of the renormalization group equation for 
the two-point vertex function and the discussion of the resulting flow equations (including their scaling behavior) 
for the coupling parameters in the framework of an explicit one-loop theory. On this basis, we shall identify the 
asymptotic critical indices for both isotropic and directed percolation, and an additional crossover exponent A. 
Finally, we shall calculate effective critical exponents for the different length scales and the wavevector dependence 
of the pair correlation function, which allows us to determine the relevant crossover scales. In the Appendices, we 
list some technical details concerning the one-loop perturbation theory, and some properties of those integrals that 
emerge after the application of Feynman's parametrization, and enter the explicit expressions for the renormalization 
constants and flow equations. 
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II. MAPPING TO A FIELD-THEORETICAL MODEL 



Following the considerations by Cardy and Sugar Q, and Benzoni and Cardy ||, we argue that our general 
percolation problem can be mapped onto a field-theoretical model. For the reader's convenience, a very brief sketch 
of the basic ideas entering the derivation of the probability measure is presented here; more details may be found in 
the literature (§,||J|]. As stated above, the central quantity of interest is the pair-correlation function G°(r2,ri) (we 
shall henceforth denote unrenormalized quantities with a superscript "0"), or probability that sites ri = (xi,ti) and 
r 2 = (x2,*2) belong to the same percolation cluster. In a first step, G°(r 2 ,ri) is represented by a sum of all such 
graphs defined on the percolating lattice which are constructed by the following rules (|J4||| (see Fig. 2): 

(1) Place oriented bonds on the lattice, in such a fashion that from each site in the diagram it is possible to reach r 2 
by following the arrows forward, and ri by following the arrows backward. (2) Closed loops of arrows are not allowed. 
(3) Insert a factor p for each bond. (4) In order to avoid multiple countings, insert a factor — i n for each vertex where 
n bonds meet. 

The second step is a more formal expression for G°(r2,ri), to be obtained by introducing commuting "ladder" 
operators a(rj) and a(rj) at each site r i: and an operation Tr, with the properties 



i a 



a 2 



(2.1) 



Tra = Tra = , Tr(aa) — 1 



(2.2) 



Denoting by P an operator projecting out the graphs with closed loops, the diagrammatic rules above can now all be 
summarized in the formula 



Introducing 



G°(r 2 ,r 1 )=Tr[Pa(r 1 ) I i [1 + pa(ri)a(rj))a(r 2 ) 

{Hnks);i,j 



A = -ln(l-p) 



the "transition probability" in Eq. (2.3) may be exponentiated according to 



]~[ [1 + pa(ri)a(rj)] = exp [Aa(r i )a(r J )] = cxp 

{links) {links);i,j 



^a(ri)Vya( r i) 



(2.3) 



(2.4) 



(2.5) 



in the final expression here, Vij — v{vi — rj) is a matrix which stems from the nearest-neighbor interaction A and 
depends on the details of the lattice structure, and v(r) is a short-ranged function. The next step in our derivation is 



a Gaussian transformation, i.e., a representation of Eq. (2.3) by a functional integral over auxiliary fields 0o and 0o 



G°(r 2 ,n) =Tr|po(n) J V[<f> ] J V[4> ] exp ( - ^ M^V^M^) 

x exp ^2 [a( r i)<Mri) + a(r.i)</>o(ri) J a(r 2 ) j . 
Now we take the continuum limit, and expand V^ 1 with respect to gradients of r; to lowest order, one has 



(2.6) 



il -^ + ^m-^dT 2 



(2.7) 
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where n c is the coordination number of the lattice. Note that we have explicitly taken account of the fact that our 
problem displays an inversion symmetry with respect to x; in the case of isotropic percolation the term oc r\ vanishes, 

-t, and hence r 



and (2.7) is even invariant under transformations t 
"time" derivative becomes irrelevant. 



-r, while in the self-affine region the second 



At last, we have to perform the operation Tr in Eq. (2.6); this may be achieved by expanding the exponential with 
respect to powers of a(r^) and a(r,). For the detailed calculations, we refer the interested reader to the Appendix 1 
of Ref . S . The final result for the pair correlation function is a sum over (m + n) point correlation functions 



°° (_ -\m+n-2 

G°(x 2 ,t 2 ;x 1 ,i 1 )= - 



?.!«.! 



Gmn( X 2^2;Xi,£i) 



(2.8) 



m,n— 1 



here the G^ n (x2, ti\ xi, £1) are defined via 

G^„(x 2 , i 2 ; x 1; h)= <0o (x 2 , i 2 ) m 0o(xi, h) n ) 

= P J V[fo] J P[0 O ] o (x 2 ,i 2 r0o(xi,ii)"exp [j[4>oM 

where the probability measure explicitly reads 



(2.9) 



J[<Po,4>o] 



cTx I dt< 0o(x,t) 



i d 2 



25o d 



dt 2 c dt 



0o (x, t) 



Uq 

' 2 



0o (x, t) 2 0o(x, t) - 00 (x, t)0o(x, 



(2.10) 



When writing down Eq. ( |2.10| ) , we have omitted nonlinear terms of higher than third order with respect to the fields 
4>q and 4>o for universality (we are interested in the scaling behavior near the percolation threshold p c ) and renormal- 
izability reasons; the neglected contributions constitute irrelevant perturbations in the sense of the renormalization 
group. Furthermore the expansion parameters of Eq. (2/7) have been renamed for convenience. For the percolation 
threshold itself one finds 



Pc 



1 



-l/n c 



and for p — > p c one has 



r <xp-p c 



(2.11) 



(2.12) 



In Eq. (2.9), the operator P now assures "causality", and projects out "acausal" diagrams, e.g.: closed Hartree loops 
containing the "response" propagator (compare Ref. 
On the basis of Eqs. (pTi 



2.10) one is now in a position to construct a perturbation expansion for the pair 



correlation function G°(r 2 ,r 1 ). In this paper, however, we shall refine ourselves to the study of the renormalization 
group equation for the two-point vertex function 

1 



(2.13) 



which already displays the correct scaling behavior and will permit the identification of the relevant critical exponents. 
Here we have defined the (d — D + 1) dimensional Fourier transformation according to 



0o(x,t) 



0o(q,^)e 



i(qx-wt) 



(2.14) 



where we have introduced the convenient abbreviation 
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(27r)° 



d D q I dui... . (2.15) 



For the renormalization of the nonlinear coupling uq, we shall also have to consider the two three-point vertex 
functions 

r o / n ,,.g a.g u\_ ^i(-f,-f;-f,-f;g,^) (9m) 

il2 V q ' W '2'2'2'2J _ G?i(q,w)G?i(f,f) 2 ' ( ' 

r0 /q^q^ \_ G; 2 (g,^;-f,-f;-f,-f) 

121 W 2' 2' 2' q ' GSiCq.w^xd,!) 2 ' 1 j 

The practical advantage of using these vertex functions instead of the correlation functions themselves is their corre- 
spondence to the one-particle irreducible diagrams within the graphical representation in terms of Feynman diagrams 
(see e.g. Ref. [0). (We remark that in the definitions of the vertex functions (2.13), (2.16), and (2.17), the Dirac S 



functions stemming from translational symmetry in x and t have been split off.) 

Returning to the "dynamical" functional pl| ( 2.10| ), it is easily seen that setting the parameter go to zero leads 



to the field-theoretical model of Benzoni and Cardy for the special case of isotropic bond percolation in d = D + 1 
dimensions j8|. Of course, cq may then be assumed to take the value 1. On the other hand, in the limit go — > oo 
and Co — > oo such that go/co remains finite, one obtains Cardy and Sugar's Reggeon field theory for the problem of 
directed percolation with t denoting the preferred direction 0] . Any finite value of go hence corresponds to a biased 
percolation problem, with go characterizing the strength of the inherent anisotropy. 

We now proceed with a simple dimensional analysis to determine the upper critical dimensions of our model in the 
different limiting cases. If we define 

[x] = [t}= A" 1 , (2.18) 

where A is a cutoff wave vector which defines the microscopic length scale of the problem, then we find for the primitive 
dimension of the stochastic fields (using [J] = A ) 

[<Mx,*)] = [<Mx,i)]=A( 2 ^)/ 2 . (2.19) 

Hence the coupling parameters acquire the following "naive" dimensions 

[r ]=A 2 , [c ]=A°, [ 5o ]=A 1 , and [u Q ] = A^' 6 ^ 2 . (2.20) 

The upper critical dimension may be identified by noting that the relevant nonlinear coupling has zero primitive 
dimension at d = d c . In the isotropic case {go — 0), the expansion parameter of the perturbation series turns out 
to be UqCo, and hence the upper critical dimension for isotropic percolation is found to be d\ = 6. In the extremely 
anisotropic limit, go — > oo (and go/co = const.), on the other hand, a simple rescaling of the fields shows that now 



the effective coupling is it 2 co/go, an d using Eq. ( 2.20 ) demonstrates that d® = 5 for directed percolation. Hence 
accompanying the crossover from self-similar to sclf-affine scaling, there is a change of the upper critical dimension. 
At first sight, this renders this crossover problem rather cumbersome, at least if one wants to use an (e = d c — d) 
expansion near d c . We shall see, however, that a procedure similar to Amit and Goldschmidt's treatment of bicritical 
points jl6| will enable us to give a consistent mathematical description of this crossover. Of course, we shall have to 
refrain from any e expansion. We would like to remark that different modifications of this (somewhat misleadingly) 
so-called "generalized minimal subtraction scheme" have been successfully employed by Lawrie [ fL8f and the present 
authors |Q9 Bl| for further interesting crossover scenarios. 
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III. PERTURBATION THEORY AND RENORMALIZATION 



On the basis of Eq. (2.10), one may derive the perturbation series with respect to uq following the common procedure 



(see e.g. Ref. |17| ). From the bilinear part of the "dynamic" functional J[4>o, </>o], one easily derives the free propagator 



G? 1(0) (q^)- ro + g2+w2/ 1 c 2_ 2icjVco , (3-1) 



and the vertices may be read off from the anharmonic part of ( 2.10] ). In Fig. 3, we depict these elements for 



constructing the Feynman graphs of our field-theoretical representation of the biased percolation problem. 

For ro — > 0, the ensuing perturbation theory is of course infrared-divergent, leading to non-trivial critical exponents. 
These anomalous dimensions are derived via studying the ultraviolet singularities of the field theory, which appear 
at the upper critical dimension d c , when the momentum cutoff A is pushed to infinity (see e.g. Ref. fl7|| ). Finite 
values are then assigned to these UV-divergent integrals through the application of a regularization prescription. We 
shall choose the dimensional regularization scheme as introduced by t'Hooft and Veltman j2^| ; here the (A — > oo) 
singularities appear as poles oc 1/ (d c — d) . 

The ultraviolet divergences may then be collected in renormalization constants and absorbed into the definition of 
multiplicativcly rcnormalizcd quantities. Thus we define the renormalized fields 



and the renormalized parameters 



t = Zl /2 fo , 4> = zl /2 4>o , (3.2) 

r=Z- 1 Z I ,(r -r 0c ) M - 2 , (3.3) 

c 2 =Z^Z- x cl , (3.4) 

g=Z- l/2 Z- l ' 2 Z g9al i- 1 , (3.5) 

u=Z-^ 2 Z u u,B 1 J 2 ^/ 2 . (3.6) 

Note that both fields are renormalized with the same Z factor, implying that Tn = Z^T^i, Ti2 = Z^^T^, and 

I 2: z,: 2 \ :i : . 



In (3J3) we have taken into account the fact that the fluctuations will also shift the percolation threshold. Further- 
more we have rendered the renormalized quantities dimensionless by introducing the explicit arbitrary length scale 
l//i, and have finally included the geometric factor 

_ r(4-d/2) 

Bd (Any/? (3J) 



in the definition of the renormalized coupling (3.£) 



In the case of crossover phenomenona, however, there is (at least) one additional relevant length scale besides the 
correlation length, given by an anisotropy or "mass" parameter describing the variation from one scaling region to 
the other. This implies the technical difficulty that both the UV and IR singularities will differ in the two distinct 
scaling regimes. In the "traditional" approach to crossover problems, one would compute the critical exponents in the 
vicinity of one of the stable fixed points; all the crossover features would then be contained in the accompanying scaling 
function as corrections to this scaling behavior. However, in general a calculation to high order in the perturbation 
expansion would be required in order to achieve a satisfactory description of the entire crossover region. Of course, 
using an (e = d c — d) expansion with respect to either of the fixed points renders the other one completely inaccessible, 
if their upper critical dimensions do not coincide. Therefore Amit and Goldschmidt's idea to incorporate the crossover 
features already in the exponent functions has proven much more successful than treating the problem on the basis 



G 



of scaling functions. The essential prescription one has to bear in mind is that the renormalization constants are not 
solely functions of the anharmonic coupling, but necessarily also of the additional "mass" or anisotropy parameter 
describing the interplay between the two different scaling regimes . 

In our case this second length scale is related to the anisotropy parameter go . For a consistent treatment of the 
entire crossover region, one thus has to assure that the UV singularities are absorbed into the Z factors for any 
arbitrary value of including go — > oo. This is not a trivial prescription, as usually the l/(d c — d) poles will be 
altered in the different scaling regimes. For the situation that we have in mind, even the value of the upper critical 
dimension is bound to change as the crossover takes place, in contrast to previously studied cases |l8|-|2l|]. However, 
we shall demonstrate that with the above stated so-called "generalized subtraction scheme" this change in the upper 
critical dimension may be incorporated into the usual formalism without any drastic changes, if one refrains from 
any e expansion about d c . The perturbation series is then an expansion with respect to the effective coupling v to 
be introduced later, which is not an a-priori small parameter, and the perturbation expansion is uncontrolled in this 
sense . If higher orders of the perturbation expansion were known, one could substantially refine the theory by a 
Borel resummation procedure. For a more detailed discussion of the question in which cases one may dispense with 
a (d c — d) expansion, we refer to work of Schloms and Dohm p4| . 

Another (minor) prize to be paid is that for the flow equations and related quantities only numerical solutions are 
accessible, and merely the limiting cases of go = and go — ► oo, respectively, allow for an analytical investigation. 
We remark that the somewhat misleading term "generalized subtraction scheme" stems from the fact that in the 
framework of an e expansion this corresponds to adding logarithms of go, which are finite in the limit e — > 0, to the Z 
factors; see the original work by Amit and Goldschmidt |T(| . It should be emphasized that the procedure described 
here is obviously applicable to a great variety of crossover problems (for some examples see Sec. V.). 

After these general statements, let us return to our explicit calculations. In Fig. 4 we have depicted the one-loop 
diagrams for T^, r" 2 , and T^. From d q 2Tn(q, 0) | g= o one infers directly the field renormalization Z^, while Z r , Z c , 
and Z g can be calculated by investigation of rn(0, 0), 9 w 2rn(0,w) | w= o, and <9 w rii(0, w) | w= o, respectively. Finally 
Z u is to be obtained from one of the renormalized three-point vertex functions at vanishing external momenta and 
frequencies (see Appendix A). All these vertex functions are investigated at finite "mass" ro = fJ 2 , in order to avoid 
complications stemming from additional infrared singularities. Thus the renormalization scale /i comes into play. 

The explicit one-loop results for the renormalization factors — being functions of the anisotropy scale go ! — read 



U0 8(6 d ~d) [ f i5(W/0-4(WM)] , (3-8) 

^^4(W.) . 0,) 
Ul 8(f-d) 6 [^(go/M)-4(go/M)] 

ulcglB^s [2liMfl) _ 4(go/M)] ( (310) 

U2 °l(6-Q 6 [4(WM)-I 3 d 3 (WM)] , (3-11) 

Z u = l- ulc f^ lUgo/») , (3.12) 

where we have introduced the abbreviations I^^go/fi) for a certain class of Fcynman parameter integrals which 
appear in the course of the calculations (compare Appendix B). 

At this point, we may study the behavior of the Z factors in the two limiting cases of isotropic and directed 
percolation, respectively. Let us investigate Z u , for example; for go — we find with Eq. ( |6.14 ) 



Z^ = 1 



Z r 



Z g = l- 
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.90 = : 



Z„ 



1 



2u 2 c B dtJl d - 6 
6-d 



while in the opposite case go — > oo the pole at d = 6 is cancelled and replaced by another one at d = 5 

ugcoS d B(l/2, (7 -rf)/2)/i rf " 5 



5o -> oo 



Z„ = 1 - 



go (5 - d) 



(3.13) 



(3.14) 



[here Eq. (6.15) has been used]. Hence our prescription provides the relevant effective coupling constants and the 
correct pole structure in both limits, and interpolates smoothly in between. 

Finally the fluctuation-induced shift of the percolation threshold results as the solution of the implicit equation 

n 2/ (6-d) 



roc 



_(d-4)(6-d) 

Note that r$ c is a non-analytical function of uq for both limits go — and go —> oo (see e.g. Rcf. PH). 



(3.15) 



IV. RENORMALIZATION GROUP AND FLOW EQUATIONS 
A. Scaling Behavior and Critical Exponents 



The renormalization group equation serves to connect the asymptotic theory, where the infrared singularities man- 
ifest themselves, with a region in parameter space, where the coupling u is finite (but not necessarily small) and 
ordinary "naive" perturbation expansion becomes applicable. It explicitly takes advantage of the scale invariance of 
the system near a critical point (i.e., the percolation threshold in our case). More precisely, we observe that the bare 
two-point vertex function is of course independent of the arbitrary renormalization scale fi: 



r?i(r ,c ,go,wo,q,^) = . 



Introducing Wilson's flow functions 



r 9 



hxZ, 



Cr = M 



d_ 

dfi 



hi : 

o r Q - r 0c 



-2 - g, + m 



dfi 



In Z r 



Cc = A* 
Cg = V 



d_ 

dfi 
d_ 
On 



Cu = MTp 

Ofl 



In — 

o c o 

ln^ 

o 9o 

ln^= 

n "0 



1, 1 1 
2 W 2^dfx 



o 

In Z c 



= -1 - <0 + Cc + /i 



d_ 

dfi 



InZ, 



a i 



rf-6 3. d 



hxZ,, 



we may transform Eq. (4.1) into a partial differential equation for the renormalized vertex function 

ru(/i, r,c,g, u, q, u) = . 



d „ d d . d ^ d 

M^- + ( r r^- + (cc— + C ff 3^- + CuU— + i 
a/x dr dc dg du 



(4.1) 

(4.2) 
(4.3) 
(4.4) 
(4.5) 
(4.6) 

(4.7) 



The symbol |o indicates that all the derivatives are to be taken at fixed bare parameters ro, go, and uo- One should 
note that, as a consequence of the generalized renormalization scheme, all the flow functions £ are functions of u, c, 
and g. 



The renormalization group equation (4.7) is now readily solved with the method of characteristics. The charac- 



teristics a(£) of Eq. (4/7) define the running parameters and coupling constants into which these transform when 
fi — > fi{€) = fit. They are given by the solutions to first-order differential equations (a = r, c, g, u) 
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i^=ma{l) , (4.8) 
with the initial conditions r(l) = r, c(l) = c, g(l) = g, and u(l) = u, namely 

a(e) = aJ? M W . (4.9) 
Defining the dimensionless vertex function Tn according to 

^n(^,r,c,g,u,c[,uj) = fi 2 f u (V,u, -, — , , (4.10) 



the solution of (4.7) reads 



^2 



(4.H) 



Here we have introduced an effective anharmonic coupling 

v = u 2 cl( 7 {g) , (4.12) 
which acquires finite values in both limits, g — > and g — ► oo. Defining the corresponding /3 function 

v , (4.13) 





the flow of the running coupling v(l) is given by the differential equation 

■ (4.i4) 

In the flow equations above, the parameter £ may be considered as describing the effect of a scaling transformation 
upon the system. Obviously, the theory becomes scale-invariant when a fixed point v*, to be obtained as a zero of 
the (3 function, 

&(«*) = , (4.15) 

is approached. The properties of Tn in the vicinity of the fixed point will yield the correct asymptotic behavior, if 
the latter is infrared-stable, i.e.: if d(3 v /dv\ , > is satisfied, for in this case the flow of the running coupling will 
reach v* for i — > 0. 



We now turn to the investigation of (4.11) near a fixed point v*; introducing the fixed-point values of the £ functions 
C — Ca(v — v*), also called the anomalous dimensions of the parameters a, we find that Tn becomes a generalized 
homogeneous function 

r u ( Mi r,c, g , W ) = pWf n (rlC, v *, q , , cV ; 2 ' 1+c) ) . (4.16) 

Using the matching condition ^ = one arrives at the following general self-affine scaling form 

rn(M, r, c, g , u, q, v) oc g '-^f n (-J^y t 1, -g^, ^ 2(g/ ^ )2 , (1 _ A) ) , (4.17) 
where we have defined /cwr independent critical exponents according to 

V± = -Q, v x = ~, z=l + C~Cg, and zA = ~Q • (4.18) 
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r]± and u± correspond to the two independent indices familiar from the theory of static critical phenomena. The 
exponent z was introduced in analogy to a dynamical critical exponent, and is in our case related to the anisotropic 
scaling behavior pq] . Finally, A is a crossover exponent, < A < 1, describing the transition from isotropic to 
directed percolation. It stems from the fact that there appear two different scaling variables for the "frequency" u> in 
Eq. ( 4.17| ). In the asymptotic limit of directed percolation, g — > oo (accompanied by c — > oo while g/c remains finite), 
the second scaling variable vanishes, and the scaling behavior is described by the three exponents ri±, v±, and z. 
Similarly, with the choice £ = (gu / ci±) 1 ^ 1+i '< i ~ i *^ Eq. 4.16 reads 

* { T Q ' 

rn(^,r, c,g,u,q, w) oc v '^Tn [ v* , — ; — rrr, 1, 



where [E5| 



Moreover, I — r leads to 



= (2 — rjx)/z , and = 



ru(/i, r,c, 5, u,q, a;) oc r 7 f u I -r v± ,—r v ^ 

fj, Cfl 



2 2 



,.-21,11 (1-A) 



with another ("static") exponent 7, related to u± and z^|| via 

7 = v ± (2 - r]±) = i/||(2 - r?||) . 
For the sake of completeness, we finally state the scaling relations || 

2/3 + 7 = 2- a = L>^ ± + i/n = 2/35 - 7 = /3(1 + 5) 



(4.19) 



(4.20) 



(4.21) 



(4.22) 



(4.23) 



for exponents a, (3, and 5, to be defined in just the same way as for ordinary continuous phase transitions |lq ]. 
In the special case of isotropic percolation (go = 0), the scaling relations become considerably simpler. For example, 



instead of Eq. ( 4.21 ) one finds self-similar behavior according to 



T u (ji,r,c,0,u,q,w) oc r^n ( l,v*, -r ",0, -5— 



(4.24) 



with two independent critical exponents r\ = — Q and f = —!/(,*. The scaling relations reduce to 7 = i/(2 — ?7) and 
2/3 + 7 — 2 — a = dv. On leaving the self-similar scaling region, a second scaling variable comes into play, leading to the 



appearance of the crossover exponent A defined above [Eqs. ( 4.17 ), ( 4.19 ), and ( 4.21 )], and eventually to anisotropic 
scaling. 



B. One— loop Results 



To one-loop order, Wilson's functions as derived from Eqs. (4.2) - (4.6) and ( |3.8| ) - ( 3.1 2| ) read 



1 



Cr 
Cc=- 



2H 

8 

vg 2 (d- 



1(7(9). 



16 



l 17V 
rd 
, J 37 



(9) 



J'' 
1 55 



iff) 



c»=-i + - 



rf-6 



Iir(9) 1(7(9) 
Iirjg)] vg 2 (d- 



1 _ -35(9) 



13v 
~16 



1(7(9) 

, ^ik(g) 
^ifr(g) 



16 



2 iM 



iJM 

1(7(9) 1(7(9) 



(4.25) 
(4.26) 
(4.27) 
(4.28) 
(4.29) 
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[here Eqs. (6.16) and (6.17) of Appendix B have been applied]. 

In the following two special situations, we may explicitly evaluate these general relations: In the case g — ► 0, we 
find v — > 2u 2 c [using (6.14)]; thus 



.9^0 



/?„ = v ( d — 6 + -v 



with the stable fixed point of isotropic percolation (for d < 6) 



v*j=-(6-d) 



(4.30) 



(4.31) 



The fixed-point values for the C functions are Q = u*/12, £ = ~ 2 + 5u*/12, C* = 0, and Q = -1 + w*/ 12 - Hence 
we find the following isotropic critical exponents [77 = rj±(vZ) and v = ^_l(vj)] 



6-d , , „ 5(6 -d) 
and i/" 1 = 2 - v 



21 



21 



In Eq. (4.17), the term linear in uj vanishes for g — > 0, and inserting Q = yields 

z(l-A) = l , 



which is equivalent to isotropic scaling in d = D + 1 dimensions; compare (4.24). The crossover exponent 

1- (6-d)/21 

2- (6-d)/21 



(4.32) 



(4.33) 



(4.34) 



provides the power law according to which the isotropic scaling region is left in favor of the anisotropic behavior to 
be discussed below. For d > 6, the isotropic Gaussian fixed point 



is stable, with the corresponding mean-field exponents 



■qx = , v x = ^, z = 2, and A = ^ 



(4.35) 



(4.36) 



The limit of directed percolation is characterized by a diverging anisotropy scale g — > 00. Using Eq. 6.15 , the 
coupling parameter becomes v — » w 2 c_B(l/2, (7 — d)/2)/g (note the additional factor l/g), and the /3 function now 
reads 



,g -> 00 : 

Therefore, for d > 5 the directed Gaussian fixed point 



V = v [ d - 5 + -v 







(4.37) 



(4.38) 



becomes stable, which again leads to the exponents of Eq. ( (4.36| ). 

On the other hand, for dimensions d < 5 the asymptotic behavior is governed by the non-trivial anisotropic scaling 
fixed point 

- 2, 



(4.39) 



Note that the upper critical dimension has changed to dj? = 5, in contrast to the isotropic case, where d\ = 6; 
physically, this reduction is due to the suppression of fluctuations by the anisotropy, and may be formally traced back 



11 



to the appearance of the factor 1/g in the expression for the asymptotic coupling. Inserting ( 4.39| ) into the results for 
the C functions, (*L — v* /8, £* = —2 + 3v*/8, Q — v* /8, and Q = — 1 + v* /4, we find the following critical indices 



_ i , 5 — d „ 5 — d 

12 1 J± ~ 4 ' " 12 ' ~~" "~ " 6 



?7_L = , i^J 1 = 2 , z = 2 , and zA = 1 , (4.40) 



which characterize the self-affine scaling of the percolation clusters with preferred "time" direction Q . The anisotropy 
is reflected in the exponent z, and A now describes the crossover from isotropic to directed percolation near the 
anisotropic percolation fixed point wfj. We have collected the four fixed points and the corresponding values for the 
independent critical exponents in Table I. 

Thus we have demonstrated that both the self-similar and the self-affine scaling behavior are within the scope of the 
present theory, at least for dimensions d < 5; for 5 < d < 6 the model is not renormalizable in the directed limit, and 
the crossover description based on extracting the UV poles may be questionable. However, in this case the asymptotics 
of the model are simply described by the mean-field exponents corresponding to the Gaussian fixed point Wq D , with 
logarithmic corrections for d — 5, and at least the qualitative features of the crossover to this Gaussian theory are 
well reproduced by our formalism (see also Refs. E3| and At any rate, this discussion again emphasizes the fact 

that no expansion with respect to a fixed upper critical dimension can be applied consistently (from our procedure a 
kind of "effective upper critical dimension" d c (g), varying with the anisotropy scale g, might be extracted). 

For the physically interesting situation d < 5, note that a smooth interpolation between the two asymptotic cases 
is obtained at every stage of the present theory, i.e., for the renormalization constants, for the C and f3 functions, for 
the fixed-point values and the scaling functions. We shall now proceed to study the entire crossover region between 
these asymptotic regimes, which can be readily done by solving the coupled set of flow equations (4.8) with (4.25) - 
(4.29) numerically. 

We start with the analysis of the flow diagram for the effective coupling "constants" v(£) and g(£) = g(£)/[l + g(£)] 
(the latter assumes values in the interval [0, 1] only, which is more convenient than the range [0, 00] of the original 
anisotropy parameter g), shown in Fig. 5. The four fixed points, as summarized in Table I, namely those for 
isotropic percolation (I), directed percolation (D), Gaussian isotropic (GI), and Gaussian directed (GD) percolation, 
respectively, determine the topology of the (v,g) flow diagram (where we have chosen d — 3 and fi — 1 for the 
renormalization scale.) The only infrared-stable fixed point is the one for directed percolation, (v* ,g*) = (§(5 — d), 1). 
All other fixed points are unstable, but as can be inferred from Fig. 5, they are more or less attractive depending 
on the initial value for the coupling constants. The flow diagram is divided into two regions by a separatrix, which 
constitutes the renormalization-group trajectory from the fixed point v{ to v^, describing the universal crossover 
from self-similar to self-affine scaling. For initial values v < v\, there are besides the stable fixed point for directed 
percolation three unstable fixed points (I), (GI), and (GD). Starting from the Gaussian isotropic fixed point (GI), the 
renormalization-group trajectories traverse regions close to the fixed points (I) or (GD), depending on the initial values 
g(l) and v(l), before they finally reach the infrared-stable fixed point for directed percolation (D). The competition 
between all these fixed points will also become apparent in the flow of the effective exponents (e.g., for the pair 
correlation function). For initial values v > i>j there is only one relevant unstable fixed point (I). 

As a result of the flow dependence of the coupling constants also the £ functions, whose fixed-point values are 
related to the critical exponents, display crossover behavior. In Fig. 6 we exemplify this flow dependence for the 
£ functions C<p(^) an( i Cr(0) respectively, where we have chosen a fixed initial value v(l) — «j for the three-point 
coupling, and a series of initial values for the anisotropy parameter [g(l) = 10~ 2 , 10~ 3 , and 10 -4 ]. 

The most important conclusion to be drawn from the crossover behavior of the £ functions is that the crossover for 
the anomalous dimension of the stochastic fields starts at values of the flow parameter which are approximately 
one order of magnitude smaller than the corresponding values for £ r . This should then also be visible as different 
crossover locations for the distinct effective exponents of the connectivity G. 
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C. Effective Exponents for the Pair Correlation Function 



In this paragraph we consider the effects of the crossover on the most interesting physical quantity, namely the 
pair correlation function. As discussed in Sec. II, the scaling behavior of the connectivity is identical to that of the 
two-point function Gn(q,u) = rn(— q, —oj)^ 1 , which has been studied in the previous subsection. 

As a first approximation, one can use the zero-loop result for the scaling function and finds 



Tii(r,q, w) = /i 



r(£) 



H 2 £ 2 fiH 2 c(£) 2 



2 1 



"g(l) 



(4.41) 



It is quite straightforward to calculate corrections to this zero-loop result for the scaling function in a perturbation 
expansion with respect to the effective coupling constant v. However, there are serious technical problems associated 
with the non-analytical dependence of the shift of the percolation threshold p c on uq [Eq. ( 3.1 5| ) ] . This is, however, 
not a particular difficulty of our theoretical approach, but a fundamental problem for any field-theoretical calculation 
at fixed dimension below d c . Schloms and Dohm |Q have shown that for a </> 4 theory the nonvanishing mass shift can 
be incorporated in the minimal subtraction approach directly in three dimensions without recourse to the standard 
(e = i—d) expansion. But, even for this "standard model" of critical phenomena there does not exist a straightforward 
perturbation method at finite external momenta which allows for a consistent treatment of the non-analytical mass 
shift at fixed dimension. Hence, at the present stage of the theory, we can treat the C functions to arbitrary loop- 
order, but have to restrict ourselves to a mean-field (zero-loop) treatment for the scaling functions. Nevertheless, one 
expects that the results for the effective exponents obtained within this "renormalizcd mean-field theory" to provide 
a reasonably good approximation. This expectation is based on the experience that amplitude (or scaling) functions 
are usually smooth and well-behaved, and enter the results less sensitively than the exponential £ functions. In fact, 
our approach was especially designed to incorporate the complete crossover behavior into the exponential functions ( a . 
This has to be contrasted with the usual scheme, where the crossover has to be inferred from high-order calculations 
of the amplitude functions (and usually does not go far beyond the determination of corrections to scaling) . 

The most convenient way to analyze the crossover behavior of the pair correlation function is in terms of effective 
critical exponents. We consider first the case r = and q = and choose to define the effective critical exponent 
V\\ eff by [compare fl4.19| )1 



2 - >7||effM 



din r n (o,o,^) 

dhiuj 



Using the matching condition 



2/ 



wg(£) 



the effective exponent is found to be 



H 2 £ 2 c 2 (£) nlc{l) 



2-r?||effW = [2 + aW] 



= 1 



dln£ 
dlnuj 



(4.42) 



(4.43) 



(4.44) 



where the factor dln£/dlnw has to be determined from 4.45 . The effective exponent 2 — r)\\ e g(£) is shown in Fig. 
7 as a function of the flow parameter £, where for the initial value v(l) — Vj, and a series of initial values for the 
anisotropy parameter, g(l) = 10~ fc with k = 2,3,4, have been chosen. Upon using the matching condition ( 4.43 ) the 
flow parameter £ can be related to the longitudinal length scale ("frequency"). If the effective exponent is plotted 
versus the running anisotropy parameter g(£), all the curves for different initial values <?(1) collapse onto one master 
curve. The corresponding plot is shown in Fig. 8. Since g(£) oc £^s near a fixed point, the scale transformation (in 
order to reduce all curves to one master curve) depends on how close the flow is to one of the four different fixed 
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points. Note that because ( g is not a constant, no simple scaling relation can be deduced for the location of the 
crossover by simply investigating any one of the two asymptotic regimes. 



Next, in the case r = and to = one can define an effective exponent r]± e s by [see Eq. (4.17)] 



dlnr 11 (0,q,0) 

2 - »7J_eff(g) = j: • (4-45) 

ding 

With the matching condition (q//i£) 2 = 1 this reduces to 

2- V±cS (£)=2 + (4£) . (4.46) 



Finally, we consider the case u> = and q = 0. Upon defining an effective exponent [compare (4.21)] 



dlnr n (r,0,0) 

7cff(r) = dhTr (4 ' 47) 

and choosing the matching condition r(£) — 1 we find 

remarkably, Eq. (|4.48|) is valid throughout the entire crossover region, and not just near one of the fixed points, where 



it becomes identical to the scaling relation ( 4.22 ) 



The flows of the above effective exponents, and that of the "effective" dynamic exponent 

ZeffW = l + CcW-C fl W , (4-49) 

are shown in Figs. 8-11 versus the scaling variable ln[flf(£)], with fi = 1, and the initial value v(l) = vf, such 
that the universal crossover starting from the isotropic scaling region is described. The most important conclusion to 
be drawn from these results is that the characteristic anisotropy scales g(^cross), where the crossover occurs for the 
effective exponents defined above, are different ! The effective exponent r]\\ e g starts to cross over from the isotropic to 
the directed fixed point value already at ln<7(£ cross ) w —0.8, whereas the effective exponent 7 c g shows this crossover 
at ln<7(^ cross ) « —0.1, and the "dynamical" exponent z c g, as well as T)± e s, at an even larger value hig(£ cross ) w 0.8. 
Note that the remarkable change of r)\\ e g is already apparent at mean-field level, where it acquires the values and 
1 in the isotropic and directed limit, respectively. However, a description of the crossovers for 7yj_ e ff, "f c &, and z e g 
requires the ( functions as least on the one-loop level, as has been achieved here. 

We remark that a precise calculation of these crossover features has not been possible up to this present work. Of 
course, the exponents for the limits of both isotropic and directed percolation have been determined to a much higher 
accuracy than is provided in our one-loop approximation p|-^0[. In fact, the relative errors of our one-loop results, 
as compared to the values given in Table III of Ref. are approximately 0.03, 0.09, and 0.21 for the exponent 7 of 
directed percolation in d — 4, 3, and 2 dimensions, respectively, and for is± one finds correspondingly 0.10 and 0.27 for 
the three- and two-dimensional cases. In the isotropic limit at two dimensions, the relative errors are 0.21 and 0.27 
for 7 and v, respectively. (We remark that the numerical values of our one-loop results are even slightly better than 
those of an additional e expansion; yet they also improve as the upper critical dimension is approached.) However, our 
aim was rather to calculate the detailed crossover properties, and we believe that the characteristic crossover scales 
<7(^cross) should not be affected too severely by, e.g., higher orders of perturbation theory. Certainly, these results, 
i.e., the approximations used in this paper, are subject to tests by both computer simulations and experiment. Of 
course, it would be very interesting to compare our predictions concerning the crossover scales of the different effective 
exponents with the outcome of such numerical simulations and/or physical experimental setups. 

Therefore we add some remarks on the interpretation of our results, and on the number of free parameters of the 
theory. One of the advantages of studying the effective exponents, is the fact that they do not depend on non-universal 
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amplitudes, i.e., on the initial values u(V), c(l) (which can be set to 1, if one starts from the isotropic scaling region), 
and the scale fi. For convenience, we have plotted our results versus the scaling variable ln[g(^)] in Figs. 8 - 11, and 
have thus also eliminated the dependence on the formal anisotropy parameter. However, in simulations or experiments 
g(l) = g is a fixed quantity, although its correspondence to a physical anisotropy measure may be rather indirect. In 
any case, if merely the universal crossover features are to be investigated, the initial value of the effective coupling 
should be chosen as v(l) — , in order to resemble the self-similar scaling behavior, and then g(l) is the only free 
parameter of the theory. In order to compare directly with our figures, one then has to solve the flow equations ( [f .8] ). 
(4.14) with Eqs. ( 4.25| ) - (4.29), and apply the relevant matching condition, which is straightforward (unfortunately, 
the scaling behavior with respect to g(l) is complex and cannot be described by a pure power law). But once g(l) 
has been determined for one of the effective exponents, it must necessarily also yield the crossover point for any of the 
others. Obviously, a quantitative study of the effective indices n\ e ff(w) and "f e s(j- ) is most promising, while a detailed 
analysis of ri± c ff(q) and z e g(£) is probably difficult, because of their comparatively small changes on approaching the 
self-affine regime. 

Possibly clearer on first sight, however somewhat less distinct concerning quantitative features, is a discussion of 
the crossover from self-similar to self-affine scaling on the basis of contour plots, e.g., for rn(0,q, w). Applying the 
matching condition 



+ 2i 



wg(£) 



= 1 



valid at the percolation threshold p = p c , Eq. (4.41) simply becomes 



rn(0,q,w) 



(4.50) 



(4.51) 



which is a monotonic function of £ (the exponential is approximately given by t-vxM*), see Fig. 9). Therefore, 
contours of constant Tn are identical to contours of constant flow parameter £, which in turn can be easily inferred 
from Eq. ( 4.50 ). In Fig. 12, we depict typical examples of such contour plots, with values of i = 10~ 4 , 10" 
10~ 2 , respectively, from which the crossover from isotropic to anisotropic scaling becomes apparent. 



and 



V. SUMMARY AND CONCLUSION 



In this paper we have investigated the crossover from isotropic to directed percolation taking advantage of a mapping 
onto a field-theoretical representation of the connectivity . From a conceptual point of view, the main result of the 
present paper is the demonstration that an extended minimal subtraction scheme is capable of dealing with crossover 
problems associated with a change in the upper critical dimension d c in the framework of renormalization group 
theory. We have exemplified the method for a long standing problem in percolation theory. But we believe that this 
approach can be applied to a wide variety of interesting physical problems. Among those of most current interest 
are the crossover between bulk and surface physics p6j ] , between mean-field and critical behavior pl[ | , crossover from 
propagating to overdamped soft modes in critical dynamics, e.g., near structural phase transitions, and others. 

Since our primary goal was to present the formalism, we have restricted ourselves to a one-loop approximation. 
We have identified the crossover exponent A and calculated effective exponents for different length scales and the 
pair correlation function. Higher loop orders for the exponential (Q functions are in principle accessible within 
our theoretical framework. However, the calculation of the wavevector and frequency dependence of the amplitude 
functions and the incorporation of the non-analytical mass shift remains an open problem for a field theoretical 
calculation at fixed dimension pi| . 

It would also be interesting to investigate the crossover from isotropic to directed percolation by numerical simula- 
tions or physical experiments, and compare them with the theoretical results of our paper. Even though the values 
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of the critical exponents in the two asymptotic regimes are not accurate (one-loop results), we expect the crossover 
behavior of the effective exponents to be qualitatively correct. Especially the predictions concerning the different loci 
of the crossover for the "static" and "dynamic" quantities could be tested by numerical simulations. Obviously, this 
would be of considerable help for estimating the quality of the proposed approach to crossover phenomena also in 
different situations, where a numerical simulation is either very cumbersome or not feasible at all. 
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VI. APPENDIX 
A. Two— and Three— point Vertex Functions to One— loop Order 

For the sake of completeness, we list some important explicit analytical results for the two- and three-point vertex 
functions to one-loop order. The corresponding Feynman diagrams are depicted in Fig. 4. Writing the free propagator 
in the form 

G° 1(0) (-q,-u;) = T -J — T , (6.1) 

where w±(q) represents the "dispersion relation" for the "elementary excitations" of our model, 



w±(q) = ~ ic o I So T y 7*0 + 0o + q 2 J , (6-2) 

one finds for r° 1 (q, oj): 

(o)= r + q 2 + lu 2 / cl + 2g iuj / c , (6.3) 

(b)= ^ J I Gu („) (k - q/2, v - u /2)G? 1 (0) (k + q/2, v + u/2) (6.4) 

= u 2 5 f 1 (qk) + ^2/2cg+2 g0 i^/c -2 g g ,q k| 

" !1 ' y/ra+^ + Wl-^f- [-^ + - + (q/2 + k) + c + (q/2-k)][-c + I ._(q/2 + k) + - + (q/2-k)] V ' ' 



k 



X 



1 



[-<u + w + (q/2 + k) + w_(q/2-k)][-u + b;_(q/2 + k) + a,_(q/2-k)] 



In the last step, we have performed the integration over the internal frequency v via the residue theorem. 

As a special case, we list the result for vanishing external wavevector, q = 0, which enters the calculation of the 
effective exponents: 



r°! 0,^ =ro + — +2,g - + ^ / _ = — — — . 6.6 

c o c o 8 J k + g$ + k 2 r + k 2 + ^ I + <? ^ 



In order to calculate the renormalization constants, we may expand (6.5) with respect to q and u) 



UqCq f 1 



8r o Jk v/r +g 2 + k 2 (r + k 2 ) ^ 

_ "qCq f 1 

co r 16 J k ^r Q +g 2 + k 2 {r Q + k 2 ) 2 



+290 ^ h-Qf , . ' (6.7) 
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+ -5" 1 



U^Co 



1 



"o c o.9o 



1 



32 Jk y/r +g% + k 2 (r + k2)2 8 J k ^r + 3o 2 + fc 2 (r + fc 2 ) 3 , 

1 MnCn /" 1 \ 



+Wi u » co 



+ 



32 J k y/r +g 2 + k 2 (r + k 2 ) 2 64 7 fe ^(ro + <^ + fc 2 ) 3 (r + fc 2 ) 



(qk) 5 



32 A V( r o +5o +k 2 ) 3 {r + k 2 )> 



+ 



3MqCq 

64 



(qk) 5 



V(ro+.g 2 + ^) 5 (^o + fc 2 ) 



Similarly, the three-point vertex function r° 2 (— q, — ui; q/2, w/2; q/2, w/2) reads 

(c) =- Uo / 2 , 

(d) = ugcg / / G° n (0) (k - q/2, „ - o,/ 2 )^ (0) (k + q/2, v + W /2)G? 1 (0) (k, «/) 



2 A 



V-O + sg + tq/a + kja l-^ + (q/2 + k) + IJ+ (q/2-k)j[-^ + ll , + (q/2 + k) + ^_ (q/2-k)] 

X ! 

[-^/2 + W+ (q/2 + k)-u, + (k)][-^/2 + W+ (q/2 + k)- W _(k)] 

-L 1 1 X 

y/ro+gZ+in/Z-k) 2 + (q/2+k)+u,_(q/2-k)J[-u, + w_(q/2+k)+u,_(q/2-k)J 

x 1 

[-^/2 + ^_(q/2-k) + ^ + (k)][-^/2 + ^_(q/2-k) + ^_(k)] 
+ v /, + g 2 +fe 2 [-^/2 + ^ + (q/2+k)- t , + (k)] 1 [- 1 ,/2 + I ._(q/2 + k)- 1 , + (k)] X 



[-u/2 + ^ + (q/2-k) + ^_ (k)J [-^/2 + ^_ (q/2-k) + c_ (k)J 

For the determination of the renormalization constant Z u , however, we merely need 
T° 2 (0, 0; 0, 0; 0, 0) - -1^(0, 0; 0, 0; 0, 0) 

_U0 ( _ ^qCQ f 1 M§CQ I" 1 

2 \ 4 Jk ^r +g 2 + k 2 (r + k 2 ) 2 8 Jk v / (^o+.g 2 + fc 2 ) 3 (ro+fc 2 ) 



(6.8) 
(6.9) 



(6.10) 



(6.11) 



B. Feynman Parameter Integrals 

Within the generalized minimal subtraction scheme described above, the determination of the Z factors after 
performing the frequency and momentum integration using the Feynman parametrization 



1 T(r + s) f 1 x^il-x)"- 1 



f 

Jo 



dx (6.12) 



A r B s T(r)T(a) J [xA + (1 - x)B] r + 
[T(x) is Euler's Gamma function] leads to integrals of the form (with odd integers m, n) 

/•l x m/2-l rg 2 j.m/2-1 

I ini9) = ] Q ( l + xg2) ( m+ n-d)/2 dx = Q~ m J o {1 + t) (m+n-d)/2 dt > ( 6 ' 13 ) 

for which one may immediately derive asymptotic relations for g = and g — > oo, respectively [B(x,y) = 
T(x)T(y)/r(x + y) denotes Euler's Beta function], 



J£ n (0)=2/m , (6.14) 
' m n — d\ r(^)r( 2 ^) 

s "-^"oo LJ - y 2"' 2 J _ r( m+ ^ rf )" 



jjm [, 9 -^( g )] = b ( V ) = • ( 6 - 15 ) 
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For the calculation of the £ functions, the following formulas are very useful: 

^„( 3 )-.g 2 ^ +2i „(.g)=/l„ +2 ( 5 ) , (6.16) 
^lL(9o/^)=(m + n-d)^Ti +2in (g /fi) . (6.17) 
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TABLE I. Fixed point values for the effective coupling "constants" v(£) and g(l) at the isotropic (I), Gaussian isotropic 
(GI), directed (D), and Gaussian directed (GD) fixed points, and the corresponding values of the four independent critical 
exponents. The other critical indices may be inferred from scaling relations. 
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Fig.l: Isotropic (a) and directed (b) percolation clusters. 
Fig. 2: Directed bonds placed on a square percolation lattice. 

Fig. 3: Propagators and vertices of the field-theoretical representation used in the text. 
Fig.4: One-loop diagrams for (a) r$ x , (b) T^ 2 , and (c) T^. 

Fig. 5: Renormalization-group trajectories in the (v,g) parameter space for the crossover 
from isotropic to directed percolation (to one-loop order) ind = Z) + l = 3 dimensions. In 
order to map the range [0, oo] of g onto the interval [0, 1], we use the variable g — g/(l + g) 
instead of the anisotropy parameter g. 

Fig. 6: Flow dependence of ($(£) (a) and C, r {£) (b), where we have chosen a fixed initial 
value v(l) = Uj for the three-point coupling and a series of initial values for the anisotropy 
parameter [g(l) = 1(T 2 , 1CT 3 , and 10~ 4 ]. 

Fig. 7: Effective exponent 2 — g\\ e e(£) for the connectivity as a function of the flow pa- 
rameter £ for fixed initial values v(l) = i>j and a series of initial values for the anisotropy 
parameter [g(l) = 1(T 2 , 1CT 3 , and 1(T 4 ]. 

Fig. 8: Effective exponent rj\\ e s(£) for the pair correlation function as a function of the 
running anisotropy parameter g(£) for a fixed initial value of the coupling constant v(l) = vf. 



Fig. 9: Effective exponent 77_i_ e ff(^) f° r the P a h correlation function as a function of the 
running anisotropy parameter g(£) for a fixed initial value of the coupling constant v(l) = Vj. 



Fig. 10: Effective exponent 7 c fr(^) for the connectivity as a function of the running 
anisotropy parameter g(i) for a fixed initial value of the coupling constant v(l) = vf. 

Fig.ll: Effective exponent z e s(£), describing the anisotropy, as a function of the running 
anisotropy parameter g{£) for a fixed initial value of the coupling constant v(l) = uj* . 



Fig. 12: Contour plots for rn(0,q,w), as inferred from Eqs. ( 4.50 ) and ( 4.51 ), with values 
of £ = 10~ 4 , 10~ 3 , and 10~ 2 , respectively. 
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